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It has long been recognized that equations describing dc transistor net- 
works do not necessarily have unique solutions. The Eccles-J ordan (flip- 
flop) circuit is an excellent example oj one for which the dc equations may 
have more than one solution. 

Only recently, however, has a comprehensive theory concerning matters 
such as the existence and uniqueness of solutions of the dc equations of 
general transistor networks begun to take shape. This paper represents 
another contribution to the evolution of that theory. 

A key concept in the development of the recent theory is the concept of a 
"P matrix." We give a generalization of that concept, showing that one 
can specify properties possessed by certain pairs of square matrices, analo- 
gous to the properties possessed by a single P matrix. Pairs of matrices 
possessing these properties are called W pairs. Use is made of this W n pair 
concept to prove results which are more general than some of the existing 
ones. We provide an extension of much of the existing theory in such a 
manner that a broader class of dc transistor netioorks may be considered. 
In particular, the new results provide one with the ability to answer certain 
questions concerning the existence, uniqueness, boundedness, and so on, 
of solutions of the equations for any network which is comprised of tran- 
sistors, diodes, resistors, and independent sources. 

I. INTRODUCTION 

Suppose a network is constructed by connecting in an arbitrary 
manner any number of transistors, diodes, resistors, and independent 
voltage and current sources. Without loss of generality, we may consider 
the network to have the canonical form shown in Fig. 1; that is, we may 
consider the network to be a multiport containing resistors and inde- 
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Fig. 1 — Canonical form of a transistor network. 



pendent sources, with transistors and diodes connected to the ports.* 
There are some fundamental questions that one should then, hopefully, 
be able to answer. For example: Do the equations that describe this dc 
network have a unique solution? With the exception of certain uniqueness 
results for a special (but none the less important) class of transistor 
networks, all of the previous explicit results in Refs. 1, 2, and 3, which 
have shown methods for obtaining answers to such questions, have 
been concerned only with the class of transistor networks for which, 
after setting the value of each independent source to zero, there exists a 
short-circuit admittance matrix (a G matrix) to characterize the linear 



* It will become apparent that the theory can also accommodate many other struc- 
tures which are of the Fig. 1 type except that the multiport contains additional 
linear elements (such as controlled sources). We do not stress this point though, 
since in the present context such elements seem somewhat unnatural. 
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multiport of Fig. 1. It is the primary purpose of this paper to show how 
that restriction can be removed. We shall show in fact that almost all 
of the previous results are but special cases of results that follow from 
a more general theory in which the assumption of the existence of a G 
matrix for the linear multiport is unnecessary.* 

Section II concerns methods for characterizing a general multiport 
containing resistors and independent sources. In Section III, we consider 
the model for a transistor. An equation for dc transistor networks is 
then developed in Section IV and, after explaining some notation in 
Section V, we develop the Wo pair concept in Section VI. Sections VII, 
VIII and IX show how the V? pair concept provides a generalization 
of the existing results concerning dc transistor networks. Finally, we 
consider an example network in Section X. 

II. LINEAR MULTIPORT CHARACTERIZATION 

A multiport having n ports (an ?i-port) is characterized by determining 
every combination of the 2n port voltages and currents that the network 
admits (see Ref. 4). We discuss here two methods of characterizing 
multiports that contain resistors and independent sources. The first 
method makes use of the f amiliar concept of a hybrid matrix. The second 
method uses a pair of matrices in a manner that was apparently first 
suggested — for multiports containing no independent sources — by 
V. Belevitch. 5 

2.1 The Hybrid Formalism 

When the value of each independent source is set to zero, for a multi- 
port containing only resistors and independent sources, the multiport 
becomes, of course, a resistive multiport. H. C. So has proved (as a special 
case of a theorem in Ref. 6) that any resistive multiport has a hybrid 
matrix description. That is, for any resistive n-port, it is always possible 
to label the port voltage and current variables in such a way that there 



* Pragmatists might argue that in any "physical" network, there will always be 
enough "stray" resistance present which, if taken into account, will guarantee the 
existence of, say, a G matrix. It seems to this writer, however, that by taking such a 
point of view, one does not obtain an entirely satisfactory understanding of matters 
(even practical matters). To know that fundamental results do not depend (if, in 
fact, they don't) upon such fortunate occurrences as these (and for many transistor 
networks this is the case) seems to be the more satisfactory situation. Furthermore, 
it should be noted that in the analysis of a physical network, to obtain a tractable 
problem, it often behooves one to neglect the presence of unimportant elements. 
Thus, it is not necessarily true that such stray resistors will always be present in the 
model of the network which the analyst desires to consider. 
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exists an integer m, ^ m ^ n, a pair of n-vectors* 
x = & , ■•• , *« , w m+] , •• • , u n ) r , 

and a real n X n matrix #, the hybrid matrix, such that the network 
admits the port variables v k , i k as the voltage and current, respectively, 
at the /cth port, for fc = 1, • • ■ , n, if and only if the vectors x and y satisfy 

y = Hx. (l) 

Thus, a resistive multiport may always be characterized by a hybrid 
matrix. 

When independent sources whose values are nonzero are present 
in an otherwise resistive multiport, a hybrid matrix will not generally 
suffice to characterize the multiport. Clearly the vectors x = y = 
(0, 0, ■ • • , 0) T which satisfy equation (1) for any matrix H do not always 
specify an admissible combination of port variables when independent 
sources are present. One might hope, however, that a characterization 
of the type 

y = Hx + c, (2) 

where c is some constant vector (whose elements are real numbers), 
might always be possible. Indeed, we are about to show that this is the 
case. There is one problem, however, that was not present in the con- 
sideration of resistive n-ports that must first be dealt with: there are ways 
to interconnect independent sources and resistors such that the resulting 
structure doesn't make sense. That is, the independent sources might 
impose self-contradictory constraints on the network. We rule out such 
possibilities by agreeing that, when we refer to "a multiport containing 
resistors and independent sources," we always assume that the multiport 
possesses the following property: 

Assumption: The linear graph that is formed by associating an edge 
with each resistor, each independent source, and each port, has no 
cut-sets containing only current source edges for which the values of 
the current sources cause a violation of Kirchhoff 's current law. Similarly, 
no circuits of voltage source edges for which the values of the voltage 
sources cause a violation of Kirchhoff 's voltage law are present. 

This assumption in no way restricts the generality of our work. We 



* We use the superscript T to denote the transpose of a vector or a matrix. Thus, 
the vectors x and y above are both column vectors. 
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are simply ruling out multiports, like the 2-port of Fig. 2, for which the 
set of admissible port voltage and current combinations is empty. 

We have worded the Assumption so that the presence of, say, a series 
connection of two 1 -ampere current sources in an otherwise resistive 
multiport does not cause the multiport to be inadmissible. We have 
done this because no violation of Kirchhoff's laws results from such 
interconnections of resistors and sources; the network is perfectly legiti- 
mate. One should be aware, however, that if "superfluous" sources are 
present in a network, it will follow that one cannot uniquely determine 
the value of each branch voltage and current in the network. That is, 
even though one might be able to uniquely determine the value of the 
voltage across the pair of 1-ampere sources, there is no way to determine 
the value of the voltage across each individual source. Aside from such 
ambiguities, it follows (see below and the proof of Theorem 1 in Ref. 6) 
that the value of all branch voltages and currents can be uniquely deter- 
mined for a multiport satisfying the Assumption, whenever the values 
of the "independent" port variables are known. 

Theorem: Any multiport containing resistors and independent sources 
can be characterized by equation (2), where H is a hybrid matrix charac- 
terization of the corresponding resistive multiport that is obtained by setting 
all independent source values to zero, and c is a vector of real numbers. 

A proof of this theorem can be constructed by incorporating a few 
simple observations and minor modifications into the arguments used 
by So in Ref. 6. We therefore simply sketch the main ideas: First, if the 
linear graph mentioned in the Assumption contains any current source 
cut-sets, then it must be the case (because of that Assumption) that 
these sources have values such that Kirchhoff 's current law is satisfied. 
That being the case, the port behavior of the multiport will clearly be 
unaltered if a sufficient number of current sources are removed (by 
coalescing appropriate nodes) to eliminate such cut-sets. A similar 
observation applies to voltage source circuits. Therefore without any 
loss of generality, we may consider the linear graph to have no current 
source cut-sets and no voltage source circuits. Next, by Lemmas 1 and 
2 of Ref. 6, it then follows that there exists a tree* for the linear graph 
for which all voltage source edges are branches and all current source 
edges are links. At each port, one of the two port variables is then desig- 
nated as "independent," the choice depending upon whether the edge 
corresponding to that port is a branch or a link. The existence of the 

* In case the linear graph is not, connected each reference to the word tree should, 
of course, be changed to forest. 
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Fig. 2 — An inadmissible n-port. 



hybrid matrix H and the vector c for the characterization (2) then follows 
in the same manner as the existence of a hybrid matrix for a resistive 
multiport follows from So's arguments. 

2.2 Belevitch's Formalism 

For some multiports, it might be that (after setting all independent 
source values to zero) a hybrid matrix exists such that the vectors x 
and y in equation (1) satisfy x = v = (vi , • • • , v n ) T and y = i — (ii , • • • , 
i n ) T . In this case the hybrid matrix is given the special name, admittance 
matrix. Similarly, if it happens that H exists such that x = i and y = v, 
then H is called the impedance matrix. For many resistive multiports, 
neither an impedance matrix nor an admittance matrix exists. It is still 
possible, however, to characterize any n-port for which a hybrid matrix 
exists in terms of the vectors v and i. Obviously, x and y satisfy equation 
(1) if and only if v and i satisfy 



[/, ! -H T ]v = [H t | -I r ]i, 



(3) 



where the n X m matrix Hi and the n X (n — m) matrix H r are defined 
by H = [Hi \ H r ], and similarly [I t \ I r ] is the n X n identity matrix. 
The characterization (3), being equivalent to equation (1), is perfectly 
adequate for any resistive ?i-port. It is, however, but a special case of a 
more general characterization due to Belevitch, namely: 



Pv = Qi, 



(4) 



where P and Q are n X n real matrices. Belevitch's characterization 
can be used for quite a broad class of networks, including some rather 
pathological ones which require dependent sources, or gyrators and 
negative resistors to realize, and for which no hybrid characterization 
exists. For example, the one-port called a norator, for which the set of 
admissible port voltage and current combinations is the set of all pairs of 
real numbers, may be characterized by [0]v = [0]i. We should note, 
however, that if one allows the aforementioned elements to be present 
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in an n-port, then even equation (4) cannot always provide a charac- 
terization. The nullator, for example, a one-port whose only admissible 
combination of port voltage and cm-rent variables is the pair (0, 0), is 
such an ?i-port. 

When an ?i-port contains independent sources it can often be charac- 
terized by the equation 

Pv = Qi + c, (5) 

where P and Q are real n X n matrices, and c is a constant vector. Clearly, 
any n-port containing only resistors and independent sources has such a 
characterization. It is this class of n-ports which is our primary concern 
in the study of transistor networks. We note, however, that equation 
(5) is adequate for characterizing a much broader class of n-ports. 



III. NONLINEAR TRANSISTOR CHARACTERIZATION 

In Fig. 3, a commonly used large signal dc transistor model is dis- 
played. It is easily verified that the voltage and current variables denned 
in that figure obey the following relationships: 

H. - — OCf 1 _ / 2 (y 2 )J 

v x I loA _ r. + r b r b \i 
v., .fJ z . L r b r c + r b _ t s J 
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Fi|f. •">— Large signal dc transistor model. 
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Each of the parameters a, and a T may assume any value in the open 
interval (0, 1). The parameters r b , r c , and r„ , which account for lead 
resistances, are sometimes omitted by device modelers (their presence 
is sometimes accounted for by including appropriate additional resistors 
in the network to which the transistor model is connected). To accom- 
modate these various points of view we specify only, therefore, that the 
values of the parameters r b , r c , and r, be nonnegative. Thus any or all 
of them may be zero. 

Depending upon whether the transistor being modeled is a pnp or 
an npn, the graph of each of the functions U an d U has one of the general 
shapes shown in Fig. 4 (at least for values of \v\ that are "not too large"). 
Often these functions are described by an equation of the form 

f k (v) = m k [exp(n k v) - 1], (k = 1, 2), (8) 

where m k and n k are appropriately chosen constants, both being positive 
for a pnp transistor, and both negative for an npn. On the other hand, 
for example, a piecewise-linear representation is sometimes specified 
for /, and j 2 ■ 

The nature of the functions j x and U for large values of |v| depends 
upon which assumptions the modeler is willing to make, and which 
effects he is interested in considering. For large negative (in the pnp 
case) values of v, for example, the graph of f k approaches — according 
to equation (8) — the horizontal asymptote i = —m k . Thus, if the 
modeler chooses to use equation (8) to describe f k for all values of v, 
the range of j k will not be the entire real line. If, on the other hand, the 
effect of ohmic surface leakage across the p-n junction is included in 
the model, the graph of the function f k will approach asymptotically 
a straight line having a small, but positive, slope. The range of such 
a function is, obviously, the whole real line. One might also wish to 
include the effect of avalanche breakdown in the reverse-biased region. 
If this is done, the graph of /* will have a shape reminiscent of that of 
a Zener diode in the v < part of its domain. 

In the forward-biased region there are also effects, particularly apparent 
for large values of v, which the modeler may or may not wish to recognize. 
For example, there is the so-called high-level injection phenomenon 
which tends to decrease the value of the forward current and which, 
using equation (8) , is usually accounted for by a decrease in the magni- 
tude of n k for large values of v. In addition, there is the effect of the ohmic 
resistance of the crystal which tends to reduce the value of forward 
current for large values of v. 

From the point of view of the device modeler, the question of whether 
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or not to include some of the effects mentioned above is often u minor 
issue. For many networks the behavior will be essentially the same 
whether or not, say, surface leakage is accounted for in the transistor 
model. From the point of view of the network analyst, however, the 
situation is somewhat different. For example, the matter of whether 
or not the functions U map the real line onto the real line can, in some 
cases, make the difference between whether or not there exists a solution 
of the network's equations. Similarly, other results that have been 
obtained recently (presented later, beginning in Section VII) also seem 
to depend upon the graphs of the functions /* having certain special 
properties. 

It seems safe to say that no matter which "special effects" are included 
(or omitted) in the description of the transistor, the functions f k may 
at least be considered to be strictly monotone increasing mappings of 
the real line into itself. For the purpose of formulating the equations 
for transistor networks, this is the only hypothesis that we shall make. 
When additional hypotheses regarding the nature of these functions 
are needed (to obtain certain results concerning properties of these 
equations) those hypotheses will be mentioned explicitly. In each case 
it will be clear that the additional hypotheses are, in some appropriate 
sense, rather weak. 

Similar remarks can be made for the diodes that are shown in Fig. 1, 
which might also be present in transistor networks. Thus, we assume 
that each diode is described by an equation of the type i = j(v) where, 
at this point, we only assume that the function / is a strictly monotone 
increasing mapping of the real line into itself. 

IV. EQUATIONS FOR TRANSISTOR NETWORKS 

Suppose we are given a dc network consisting of transistors, diodes, 
resistors, and independent voltage and current sources, connected to- 
gether in an arbitrary manner. Let there be n transistors and d diodes. 
Clearly, there is no loss of generality if we consider the network to be of 
the type shown in Fig. 1. Using the results of Section III, we may describe 
the nonlinear devices in the network by the equations 

y = TF{x), x = £ - Ry, (9) 

where T = dh\g[T l , T 2 ], with 7\ a block diagonal matrix with n2 X 2 

diagonal blocks of the form 

1 << 

1 — a r 



(*) 



-— «/ 



(*-) i 



for k = 1, ••• ,n, (10) 
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Fig. 4 — General shape of the functions /i and/gj (a) pnp transistor, (b) npn transistor. 

and T 2 the cl X d identity matrix. Also, R = diag [Rt , R 2 ], with R x a 
block diagonal matrix with n 2 X 2 diagonal blocks of the form 



+ tf 






+ r 



k = 1 



(11) 



and /2 2 the cl X d matrix whose elements are all zeros. The function F 
has the form F(x) = [/i(xi), • • • , f2n+d(x 2n +d)] T , where each of the f k is a 
strictly monotone increasing mapping of the real line into itself. 

Using the results of Section II, the effect of the linear multiport in 
Fig. 1 is to constrain the vectors of port variables, # and y, to obey the 
relationship 



Pz = -Qy + c, 



(12) 



where P and Q are (2n + d) X (2?i + d) real matrices and c is a real 
(2n + d)-vector. The minus sign appears in equation (12) as a consequence 
of having chosen the reference direction for the port currents (the elements 
of the vector y) to be opposite to that which is usually assumed. 

By using equations (9), we may easily eliminate the variables £ and 
y from equation (12), resulting in the equation 



(PR + Q)TF(x) + Px = c. 



(13) 



The central problem in determining the values of all branch voltages 
and currents in a dc transistor network is the determination of a solution 
of equation (13) . The rest is relatively straightforward, for if x is a (unique) 
solution of equation (13), then the (unique) vectors £ and y, such that 
equations (9) and (12) are satisfied, may immediately be computed 
from equations (9). 

Since the matrix T is nonsingular, it follows that whenever either 
(PR + Q) or P is nonsingular, equation (13) can be transformed into, 
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respectively, one of the equations 

F(x) + Ax = b, (14) 

AF(x) + x = b. (15) 

The first of these equations has been studied rather extensively (see 
Refs. 1-3 and 7) and for most of the results obtained there, it can be 
shown that parallel results are possible for equation (15). Both of these 
equations, however, are but special cases of the equation 

AF(x) +Bx = c, (16) 

which accommodates equation (13) directly. It is, therefore, this equation 
to which we shall now direct our attention. It will be shown that most 
of the results which have been obtained to date for equation (14) have 
rather natural (though not obvious) extensions to equation (16). It is 
important that such extensions be possible because one is often forced 
to deal with equations like (16) in the analysis of transistor networks. 
Clearly, this is the case whenever both of the matrices (PR + Q) and 
P of equation (13) are singular— and this can easily happen (for example, 
if the matrix R contains all zeros, then it will happen whenever there 
exists no admittance matrix nor impedance matrix for the linear multi- 
port of Fig. 1). 

V. NOTATION 

The following notation shall be used throughout the remainder of 
the paper: For each positive integer n we denote by W the n-dimensional 
Euclidean space, the elements of which are ordered n-tuples of real 
numbers, which we consider to be column vectors. The origin in E" is 
denoted by 0. If z = (rc x , • • • , x n ) T and y = fa , • • • , y,) T are elements 
of W we denote their inner product by (x, y) = 2"=i x hVk • The norm of 
each x £ E n is denoted by ||x|| = (x, x) h . 

If A is an n X n matrix, then for A; = 1, • • • , n, A k denotes the fcth 
column of A. A principal submatrix of a square matrix A is any square 
submatrix of A whose main diagonal is contained in the main diagonal 
of A. A principal minor of A is the determinant of any principal sub- 
matrix of A. If D is a diagonal matrix, then D > means that each 
element of the main diagonal is a positive number; similarly, D ^ 
denotes that each element of the main diagonal is nonnegative. We 
denote the n X n identity matrix by either 7„ or, when the dimension 
is unimportant or is clear from the context, simply by I. The direct 
sum of two matrices A, Bis denoted by A © B. A square matrix of real 
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numbers A is said to be strongly row-sum dominant if its elements 
a,,- satisfy a,, > /, lali |a,,| for i = 1, • ■ • , n. 

If / is a real valued function defined on E 1 then / is said to be monotone 
increasing if for all x < y it follows that j(x) ^ /(y). We say that / is 
strictly monotone increasing if f(x) < /(y) for all x < y. For each positive 
integer n, we denote by SJ" that collection of mappings of E n onto itself 
defined by: F £ SF" if and only if there exist, for i = 1, • • • , n, strictly 
monotone increasing functions /,• mapping E x onto E 1 such that for 
each x =(*!,-•• , rr„) r G ET, F(x) = [/,(*,), ■ • • , U(x n )] T . 

VI. PAIRS OF MATRICES OF TYPE %V 

Many of the recent results referred to above, concerning equation 
(14), have relied heavily upon certain properties that a matrix is known 
to possess whenever it is a member of a class of matrices that has been 
given the name P . In a similar way the results that follow rely upon 
useful properties that are possessed by certain pairs of matrices. We 
shall define a class, the elements of which are these pairs of matrices, 
and give it the nameW . 

The class of matrices called P was defined by M. Fiedler and V. 
Ptak. 8 They proved that the following properties of a square matrix of 
real numbers, A, are equivalent: 

(i) All principal minors of A are nonnegative. 

(ii) For each vector x 5* there exists an index k such that x k ^ 

and x h (Ax) k ^ 0. 
(Hi) For each vector x ^ d there exists a diagonal matrix D z ^ 

such that (a;, D x x) > and (Ax, D x x) ^ 0. 
(iv) Every real eigenvalue of A, as well as of each principal submatrix 

of A, is nonnegative. 
Sandberg and Willson proved that another property can be added to 
this list of equivalent properties, 2,3 namely: 

(v) det (D -\- A) 7± for every diagonal matrix D > 0. 

The class of all matrices possessing one (and hence all) of the above 
properties is called P . 

We shall now state a theorem which provides a useful generalization 
of the concept of the class of P matrices. 

Definition: For each pair of n X n matrices (A, B) we shall denote by 
Q(A, B) the collection of all the n X n matrices that can be constructed 
by juxtaposing columns taken from either A or B while maintaining 
the original relative ordering of the columns. Thus, M £ G (A, B) if 
and only if for each k = 1, • • • , n, either M k = A k ov M k = B k . 
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Obviously e(A, B) contains 2 n matrices (for certain pairs (A, B) — 
namely for those having A k = B k for one or more values of k— it can 
happen that two or more matrices in Q(A, B) are identical). 

Definition: The pair of n X n matrices (M, N) is said to be a comple- 
mentary pair taken from e (A , B) if and only if both M and N are members 
of Q(A, B) and for each k = 1, • ■ ■ , n, either M k = A k and N k = B k , 
or else M k = B k and N k = A k . 

It is obvious that (A, B) is a complementary pair taken from Q(A, B). 
It is also clear that Q(A, B) = e(B,A) and, moreover, that if (M, N) is 
any complementary pair taken from Q(A, B), then Q(M, N) = Q(A,B). 
Furthermore, for each M G Q(A, B) there exists N EQ(A,B) such that 
(M, N) is a complementary pair. 

Theorem 1: The following properties of a pair of n X n matrices of real 
numbers (A, B) are equivalent: 

(i) del (AD + B) ^ for every diagonal matrix D > 0. 

(ii) There exists a matrix M E Q(A, B) such that del 1^0 and such 

that del M-det N ^ for all N G 6(A, B). 
(Hi) For each vector X9*d there exists an index k such that either (A T x) k 7± 

or (B T x) k j± 0, and such that (A T x) k (B T x) k ^ 0. 
(iv) For each vector x ?± 6 there exists a diagonal matrix D x ^ such 

that either (A T x, D x A T x) > or (B T x, D x B r x) > (that is, 

such that (A T x, D z A T x) + (B T x, D z B T x) > 0), and such that 

(A T x, D x B T x) ^ 0. 

(v) For each complementary pair of matrices (M, N) taken from 

Q(A, B), each real value of X that satisfies det(M — \N) - is 

nonnegative. 
(vi) There exists a complementary pair of matrices (M, N) taken from 

e(A,B) such that M~ l N G P ■ 
(vii) There exists a matrix M G C(A, B) such that det M ^ 0; and, 

for any complementary pair of matrices (M, N) taken from e(A, B) 

with det M ^ 0, il/ _I iV G ^o . 

In this paper, we do not make use of properties (Hi), (iv), or (w) of 
Theorem 1. The proof that the remaining four properties are equivalent 
is given in the Appendix. A complete proof of Theorem 1 is given else- 
where. 9 

Definition: The class of all pairs of matrices which possess one (and 
hence all) of the properties listed in Theorem 1 is called W . 
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To see that properties (i) and (ii) of Theorem 1 are in fact generaliza- 
tions of the previously mentioned properties (v) and (i), respectively, 
that define P is a simple matter. It happens that for any n X n matrix 
B the pair (J„ , B) E^o if and only if B G P . (This follows from pro- 
perty (pit) of Theorem 1.) With our attention restricted to pairs of 
matrices of the type (J n , B), it is clear that property (i) of Theorem 1 
is equivalent to property (v) which determines those matrices B that 
that are in P Q . Concerning property (ii) of Theorem 1, an arbitrary 
matrix N £ G(I„ , B) is either the matrix /„ or else, a matrix formed from 
B by replacing some of the columns of B by the corresponding columns 
of I n . Consequently, det N = det B N where B N is the principal sub- 
matrix of B that is formed by removing from B the columns that are 
not present in N and then removing the corresponding rows. Hence, 
since det /„ 5^ 0, we may take I n to be the matrix M in property (ii) of 
Theorem 1, and observe that this property then becomes: det B N ^ 
for all N G C(/ n , B). It is now clear that this property is equivalent 
to the property (i) that defines the class of P matrices. (Note that there 
are exactly 2 n — 1 principal minors for each n X n matrix, and that 
the set e(/„ , B)\{I n } contains exactly 2" — 1 members.) 

VII. THEOREMS ON EXISTENCE AND UNIQUENESS 

7.1 First Existence and Uniqueness Theorem 

The following theorem, which is proved in Ref. 2, provides a necessary 
and sufficient condition for the existence of a unique solution of equation 
(14) for all F that are strictly monotone increasing "diagonal" mappings 
of FT onto FT and for all b G FT. 

Theorem 2: If A is an n X n matrix of real numbers, then there exists a 
unique solution of equation (14) for each F G 5" and for each b G FT if 
and only if A ^ P„ . 

Using this theorem along with the results of Section VI we can prove 
the following (more general) theorem. 

Theorem 3: If A and B are n X n matrices of real numbers, then there 
exists a unique solution of equation (16) for each F G 5" an d each c G F" 
if and only if (A, B) G Wo • 

Proof: (if) Let (A, B) G "Wo • Then, by Theorem 1, there exists a com- 
plementary pair (M, N) taken from e(A, B) such that M~ l N G P • 
For each P= &(•), -•• ,/.(.)]' G 5" let (7= [ 9l (-), •■■ , g n (.)f denote 
the mapping (also in %") defined by 
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it. 



»(.)- «■> if ^ = A - for *-!,■ 
|fr»(.) if M k * A t , 

Clearly, the vectors x and y satisfy 

AF(x) + Bx = MG(y) + Ny 

if they satisfy the relation 

Vk = \ X " if **> = **> for /c = l, ••• ,n, (17) 

(f h (x h ) if M * ^ A k , 

and since this relation defines a homeomorphism of E n onto itself, it 
follows that there exists a unique solution of equation (16) for each 
c E E" if there exists a unique solution of the equation 

MG(y) +Ny = c (18) 

for each c G BT. But, that this is so follows immediately from Theorem 
2 and from the fact that M~ l N E P» ■ 

(only if) Suppose (A, fl) $*W„ . Then, by Theorem 1, there exists a 
diagonal matrix D > such that det(AZ) + B) = 0. Choosing F(x) = Dx, 
we have F £ IF", while equation (16) does not have, with this choice of 
F, a unique solution for all c £ W. □ 

There are corollaries to Theorem 2, given in Ref . 2, that also may be 
generalized in a similar manner. For example, the following result is a 
generalization of an important special case of Corollary 1 of Ref. 2; 
it shows that the condition (A, B) E W is still sufficient to insure the 
uniqueness of a solution of equation (16) (if a solution exists) even 
when the mapping F is not onto. 

Theorem 4: // F(x) = [/i(.r,), ■ ■ ■ , fn(x n )] T , where each j k is a strictly 
monotone increasing mapping oj E l into E 1 , and if (A, B) E "Wo , then 
there exists at most one solution oj equation (16) jor each c E E*> 

Proof: Suppose that, for some cEFTjX 1 and x 2 are solutions of equation 
(16) with x 1 - x 2 * 6. Then, A[F(x l ) - F(x 2 )] + B(x l - x 2 ) = d. But 
then, since F is a strictly monotone increasing "diagonal" mapping, 
there exists a diagonal matrix D > such that F(x l ) — F(x 2 ) = 
D(x l - x 2 ), and hence (AD + B)(x l - x 2 ) = d. Since x l - X* * 6 it 
follows that det(AD + B) - 0, which implies that (A, B) $«W . □ 

7.2 A Nonuniqueness Theorem 

From the proof of the "only if" part of Theorem 2 (given in Ref. 2) 
it follows that whenever A ^ P , there exists a mapping F E $ n and a 
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vector b £ E n such that equation (14) has more than one solution. On 
the other hand, even if A (£ P , if the mapping F £ 5" is "fixed," then 
it is easy to see that the nonuniqueness of solutions of equation (14) 
need not necessarily follow for any b £ iTftake F(x) = x and Ax = — 2x, 
for example]. I. W. Sandberg has shown, 10 however, that if one assumes 
that the "fixed" mapping F has another special property, rather than 
assuming that F £ 5", then the nonuniqueness of solutions of equation 
(14) follows, for some b £ E n , whenever A (£ P . Moreover, he has 
shown that under these hypotheses and for any 8 > 0, there exists some 
b £ E n such that equation (14) has two solutions, x and y, which satisfy 
\\x — y\\ = 8. The special property that F is assumed to have is given in 
the following definition (in words, the property is: that it be possible 
to draw a straight line having any given positive slope, and any given 
length, between some pair of points on the graph of each of the functions 
/*)• 

Definition: For each positive integer n we denote by 8" that collection 
of mappings of E n into itself defined by: F £ 8 n if and only if there exist, 
for k = 1, • • • , n, continuous functions f k mapping E 1 into E l such that 
for each x £ E n , F(x) = [fifa), • • • , f n (x n )] T , with each of the f k satisfying, 
for all /3 > 0, 

inf {/*(« + 0) - /,(«): -co <«<«>!= 0, 

sup {f k (a + |8) — f k (a): -co < a < cc } = oo . 

By using Theorem 1 it is possible to prove the following generalization 
of Sandberg's result: 

Theorem 5: Let F £ S n , let (A, B) £JE W be a pair of real n X n matrices, 
and let 8 be a positive constant. Then, for some c £ E n there exist solutions 
of equation (1G), x and y, satisfying \\x — y\\ = 8. 

Proof: Since (A, B) (£ W there exists a diagonal matrix D = 
diag(dj , • • • , d n ) > 0, such that det(AD + B) = 0. Therefore, there 
exists x* £ E", with ||.c*|| = 8, such that (AD + B)x* = 6. Since F £ 8" 
there exists x £ E n such that 

f k (x k ) - f k (x k - x%) = x* k d k , for k = 1, • • • , n. 

Let c = AF(x) + 2?.r, and let ?/ = x — x*. Then 

4[F(a>) - F(y)] + J5(.r - y) = A[F(x) - F(x - x*)] + Bx* 

= (AD + B)x* = e. n 

For a mapping F to be a member of 8", it is not necessary that F £ $". 
It follows from the above definition of 8" that F £ 8" implies that each 
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of the functions /* is a monotone increasing function from E l onto some 
interval in E l whose length is infinite; the /* need not, however, be 
strictly monotone increasing, nor onto E\ For those F G 6" for which 
each of the functions f k is strictly monotone increasing, we have the 
following corollary to the two preceding theorems. 

Corollary: LetF(x) = [fiM, • ■ ■ ,/»(«„)] T E S," and let each of the functions 
f K be strictly monotone increasing. Then there exists at most one solution 
of equation (16) for each c G #" */ and only if (A, B) G *W . 

VIII. RESULTS ON CONTINUITY AND BOUNDEDNESS 

For many systems whose behavior is described by an equation having 
the form (16), the vector c may be regarded as the system's input and 
the vector x may be regarded as the system's response or output. Those 
properties that one might expect well-behaved systems to possess are 
likely to include continuity and boundedness. Thus, one might expect (*) 
"small" changes to result in the value of the system's output when 
"small" changes are made in the value of the system's input, and (u) a 
bounded sequence of input vectors to yield a bounded sequence of 
outputs. We now show that such properties are indeed possessed by 
the type of system that is the main concern of this paper. 

8.1 Continuity 

When the n X n matrix A is a member of the class P„ and the mapping 
F G 5", it follows that the solution x of equation (14) is a continuous 
function of the (input) vector b. a Using this fact, it is easy to prove the 
following theorem. 

Theorem 6: For each F G ff* and each pair of n X n matrices (A, B) G W 
the solution x of equation (16) is a continuous function of the vector c. 

Proof: Proceeding as in the "if" part of the proof of Theorem 3, we see 
that the theorem follows immediately from the facts that equation 
(17) is a homeomorphism and that the aforementioned result guarantees 
that y, the solution of equation (IS) , is a continuous function of c. □ 

8.2 Boundedness 

In Ref. 2 a theorem (Theorem 5) is proved which shows that, when 
F G $" and A ^ P , bounds can be obtained for the solution of equation 
(14) whenever bounds for b G E" are given. The proof of a more general 
theorem concerning equation (16) can be constructed quite easily by 
using that theorem, and by using the same technique that was used in 
the proof of the preceding theorem, along with the trivial observations: 
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(?') For any nonsingular n X n matrix of real numbers, M, and any 
real numbers a, ^ /3, ■ , i = 1, • • • , n, there exist real numbers, 
«' ^ &'i , i = 1 , • ■ • , n, such that when each of the components 
d of the vector c satisfies a, ^ c, ^ /3, , it follows that 
a', ^ (M~ l c)i g j8J , for i = 1, ■ • • , rc. 

(«) For any given real numbers 7. ^ 5, , i = 1, • • • , ?i, there exist 
for the homeomorphism (17), real numbers 7 J ^ 5' , i = 1, • • • , n, 
such that whenever x, y satisfy equation (17) with 7, ^ ?/, ^ 5, , 
for i = 1, • • • , n, it follows that 7' ^ a;,- ^ 5J , for i = 1, • • • , n. 

The more general theorem, whose quite obvious proof is omitted, 
is the following: 

Theorem 7: Let F £ 5", let (A, B) GW be a pair of n X n matrices, and, 
for i = 1, • • • , n, let a, ^ /3, fee given. There exist, fori= 1, ■ • • , n, real 
n umbers 7, ^ 5, such that for any c = (ci , • • • , c,,) 7 ^ £ E n withai ^ c,- ^ /3, 
/or * = 1, ■ ■ ■ , n, if x satisfies equation (16), //ie?i 7, ^ a;, ^ 5, for i = 
1, • •• ,n. 

According to Theorem 7, (A, B) £ %V is a sufficient condition for a 
bounded sequence of vectors c to yield a bounded sequence of solution 
vectors of equation (16), for all F £ ff". The following theorem shows 
that (A, B) G"Wo is also a necessary condition. 

Theorem 8 : If (A,B) is a -pair of real n Xn matrices, then (A,B) £ Wo if 
and only if 'for each F £{F n and each unbounded sequence of points x l ,x 2 ,x 3 , • • ■ 
in FT, the corresponding sequence c 1 , c 2 , c 3 , • • • [c = AF(x k ) -\- Bx k , 
k = 1, 2, 3, • • ■] is unbounded. 

This theorem, which is a generalization of Theorem 4 of Ref . 2, can be 
proved in a manner which is a quite obvious generalization of the proof, 
given there, of that theorem. Thus, an appeal to Theorem 7 proves the 
"only if" part, and the "if" part is proved by assuming that (A, B) ^W 
and then choosing the same land of mapping F £ JF" as was chosen in 
Ref. 2, for which an unbounded sequence of vectors x k yields a bounded 
sequence of vectors c. 

DC. COMPUTATION OF THE SOLUTION 

A. Gersho 7 has shown that whenever F £ 5" C\ C l (that is, whenever 
each of the functions /* is a continuously differentiate strictly monotone 
increasing mapping of the real line onto itself), it is possible to compute 
the solution of equation (14), for any A £ P and any b £ FT, by making 
use of a gradient descent algorithm due to A. A. Goldstein. n The following 
theorem extends this result to the class of equations of the type (16). 
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Theorem 9: Let M be an arbitrary positive definite symmetric matrix, 
and letQ : E" — > E 1 be defined by 

Q(x) = [AF(x) + Bx - c] T M[AF(x) + Bx - c], 

where F E S" H C\ (A, B) E W , and c E FT. For each x E E" and each 
y ^ let 

[Q(x) - Q[x - yVQjx)] 



9(z, y) = 



7 > 0; 



7||VQ(.r)H 2 

1, 7 = 0; 

where VQ(z) denotes the gradient of Q at the point x. Then, if 8 is any real 
number satisfying < 8 ^ \, and if x° is an arbitrary point in E n , the 
sequence {&': h = 0, 1, 2, • • • } converges to the solution of equation (16), 
where {for k = 0, 1, 2, ■ • •) the x k satisfy 

z* +1 = x * - y k VQ(x k ), 

each 7* being any real number that satisfies 8 = g(x , y ) ^ 1 — 8 if 
g(x k , 1) < 8,ory k = lifg(x k , 1) = 5. 

Proof: This proof uses generalizations of some of the ideas in Ref . 7 and 
relies ultimately upon the Goldstein algorithm. 11 

We first remark that the sequence {x k } is well-defined: It is easy to 
show (see the first part of the proof of Theorem 1, p. 31, Ref. 11) that 
for each x £ FT, g(x, •) is a continuous function on [0, °o). This being 
the case, it is clear that if g(x k , 1) < 8, then for each £ in the interval 
[S, 1) — and, in particular, for each £ in the interval [5, 1 — 5] — there is 
some 7* in the interval (0, 1) such that g(x", y k ) = £. 

Let 5 = \x E FT: Q(x) = Q(x°) } . Using the fact that M is a positive 
definite symmetric matrix, and using the fact that F E ff", (A, B) £*Wo 
implies that || AF(x) + Bx || -* «> if and only if || x || -> «> (Theorem 8) 
we have that the set S C E n is bounded. By continuity of Q, S is closed. 
Thus, S is compact and, therefore, the gradient VQ (which is continuous 
on E n , since F E C 1 ) is unif ormly continuous on S, and VQ is bounded 
on S. Also, Q is bounded below on S. [Indeed, we have Q = on E" and 
by the existence and uniqueness theorem, Theorem 3, there exists 
exactly one point x* (x* E S) at which Q(x*) = 0.] 

It is easily verified that, for each x E E n , 

VQ(s) = 2(AD X + B) T M[AF(x) + Bx - c], 

where, for k = 1, • • • , n, the /cth diagonal element of the diagonal matrix 
D x > has the value of the derivative of the function /* , evaluated at 
the point x h . Since (A, B) E Wo implies that det (AD X + B) ^ 0, and 
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since det M ^ 0, it follows that VQ(.r) = 6 if and only if x is the solution 
of equation (16). 

In view of the above, it follows directly from Goldstein's theorem that 
the sequence \x*] converges to the solution of equation (16). □ 

Other methods of computing the solution of equation (16), in certain 
cases, also exist. If one performs a transformation of the type (17) on 
the independent variable x (in theory this can always be done) then the 
solution of equation (16) can easily be computed by first computing the 
solution of an equation of the type G(y) + M~ 1 Ny = il/ _1 c, where 
G £ JF" and M~ 1 N £ Po ■ Methods of computing the solution of certain 
equations of this type may be found in Refs. 1-3. 



X. EXAMPLE 

With the aid of the modern computing facilities that are commonly 
available today, it is clearly a rather routine matter to obtain an equation 
of the type (16) for any given transistor network. Moreover, it is not 
unfeasible, even for networks of moderately large size (say, up to 4 or 
5 transistors), to consider the straightforward evaluation of the 2" 
determinants specified in property (2) of Theorem 1, and thereby resolve 
the issue of whether or not the matrices involved in the equation are 
a Wo pair. Due regard would of course have to be paid to the matter of 
performing sufficiently accurate computations. 

On the other hand, even without the aid of a computer, it should often 
be possible to use a little ingenuity and a few devices* to reduce the 
computations involved in the application of the above theory to many 
specific problems to a point where they will just about fit onto the back 
of an envelope. Consider, for example, the following analysis of a three- 
transistor network: 

For the network of Fig. 5, the voltage and current variables defined 
there must satisfy the following equations: 



(19) 



*1 






Ufa) 


/, 




Ufa) 


h 




r (l) i 
i _ i 


Ufa) 


h 




1 y,« ] 
1 _ 1 


Ufa) 


h 




1 1 rpl3) 


Ufa) 


U 






.Ufa) 



* According to It. Bellman: "a device is a trick that works at least twice." 12 
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Fig. 5 — Example of a three transistor network. 
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(20) 



l-i-j l -1 I G 2 G 2J 

where (we are using the transistor model of Fig. 3, with r b = r„ = r, = 0) 
each of the 2 X 2 matrices T {k) , k = 1, 2, 3, is of the form (10). A hybrid 
characterization has been used for the linear part of the network. As 
indicated in equation (3), this hybrid characterization can easily be 
converted into a characterization of the Belevitch type. Thus, denoting 
the 3 X 3 blocks of the hybrid matrix in equation (20) by H u , H l2 , H 2 i , 
H 22 , in the usual manner, one obtains 



/ -Hv. 
-H 2 ., 



H 21 





-I 



i + c, 



(21) 



where v — (v i , v.< , v 6 , v 2 , v 4 , v 5 ) T and * is similarly defined. We could now 
simply reorder the columns of each matrix in equation (21) in such a way 
that the resulting equation would have the same form, except that the 
subscripts on the components of the vectors v and i would occur in the 
natural order (1, 2, 3, 4, 5, 6) and then use that equation, along with 
equation (19), to produce an equation of the type (16) for our network. 
In this example, though, it's probably easier to reorder the rows and 
columns of the matrix T (recall, T = T (1) © T (2) © T (3) ) to obtain 
from equation (19) an equation that is compatible with equation (21). 
Thus, 

r t -pi 

F f(v), (22) 

-Q I J 



i = 



1734 
where 

and 
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F(») = [/.(»>), 1M, UM, U(v 2 ), 1M), Ufa)] 7 ", 



P = diag [««>, a™, a? '], Q - diag [a} 1 ', «< 2 \ «<«]. 
Eliminating ?' from equations (21) and (22), we obtain 



LF 21 -/ 



I -P 



TO 



/ -H 12 
-77 22J 



y = c. 



(23) 



Note that since detH n = detH 2 2 = 0, it is impossible to put this equation 
into either of the forms (14) or (15). Clearly this would be the same situa- 
tion no matter which ordering of subscripts was chosen for the components 
of v. The cause of the difficulty is simply the fact that neither an impedance 
matrix nor an admittance matrix exists for the linear part of our network. 
Let us determine whether or not the pair of matrices 



Hu 



H 21 -/__- 



I -P 



I 



-# 22J 



is a %*?„ pair. We shall try to verify property (1) of Theorem 1. Let 5, , • • • , 5 6 
denote arbitrary positive real numbers, and let A T = diag (5. , 5 2 , 8 3 ), 
A IT = diag (5 4 , S B , S 6 ). We wish to show that 



det 



\\H U 

IA. -i 



/ -P A~r 







+ 



/ -H lt 



_-Q 7 JL A„j L0 -H 22 _ 



^ 0. 



By multiplying the above matrix on the left by the (nonsingular) matrix 
diag (7 : , , — 7 3 ) and then multiplying on the right by diag (A, , 7 3 ), we 
obtain the equivalent statement: 



det 



-H 2 7- Q 



—H 12 — H U P A u 
H 22 + (7 + H n P) A, J 



^ 0. 



The 3X3 submatrix in the upper left corner is nonsingular and diagonal. 
The 3X3 submatrix in the lower left corner can be diagonalized by 
performing a single elementary row operation on the matrix; namely, by 
subtracting l/(5 2 + Ri) times the second row from the fourth row. 
Having done this, our problem reduces to one of showing that 
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It is easy to verily that whenever det A n 9^ 0, then 

in in 



det 







*o. 



i 21 i 22 . 

if and only if det (A n — Ai\A~^A x2 ) ^ 0. In our case both A u and A 2l 
are diagonal and hence we can immediately reduce our problem to: 

g£ « j+(1 _ 0as+ i^ (1+a :" ff , 3j) Gi 

G, G 3 +(l-a; J, )«.+^^(l+a;"R,«.) 

It is obvious that the above determinant is always positive. First, note 
that every term in the matrix is nonnegative except, possibly, the 
(1, 2) term, which may be either positive or negative (or zero). In the 
event that the (1, 2) term is positive (or zero), we have l/(5 2 + R t ) ^ 
(1 — a ( r 2) 5 2 <5 5 ) / (<5 2 + Ri), and hence we observe that the matrix is 
strongly row-sum dominant. This implies that its determinant is positive. 

In the event that the (1, 2) term is negative, we do not necessarily 
have dominance; however, considering an expansion of the determinant 
along its first row we sec that, because of the assumption that the (1, 2) 
term is negative, the value of the determinant is computed as the sum of 
two positive terms. 

We have thus shown that, no matter which (positive) values are 
assigned to R t , R 2 , R 3 , or which values the transistor's current gains 
assume [0 < a ( f k) < 1,0 < a l r k) < 1], the pair of 6 X 6 matrices that appear 
in equation (23) is a W pair. Thus, all of the results concerning a solu- 
tion's existence, uniqueness, continuity, boundedness, and so on, hold 
for this equation. 

XI. ACKNOWLEDGMENT 

The author has benefited from discussions with his colleagues I. W. 
Sandberg and H. C. So. 



1736 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1970 

APPENDIX 

Proof of Part of Theorem 1 

In this appendix we prove the equivalence of properties (i), (ii), (vi), 
and (vii) of Theorem 1, which define the class of pairs of matrices W . 
We omit the proof of the equivalence of the three remaining properties, 
since those properties are not referred to in this paper. A complete proof 
of Theorem 1 is given elsewhere. 9 We begin by proving a useful lemma: 

Lemma 1: For each positive integer n the polynomial 

(c ) d, d 2 ■■■ d n + (cj d.ch ■■■ dn-j + • • • + (c„) d 2 • • • d n 

+ (c n+1 ) dt d 2 ■■■ d n - 2 + • • • + (c n(n+1)/2 ) d 3 ■■■ d n + ■■■ + (c 2 „_,) 

in the n variables d t , d 2 , •■■ , d n is nonzero for all positive values of the 
variables if and only if at least one of the coefficients c , • • • , c 8 ._i is nonzero, 
and all nonzero coefficients have the same sign. 

Proof: (By induction) For n = 1 the statement is obviously true. 
Let N be a positive integer. Then any polynomial of the above type in 
N + 1 variables, (c )di • ■ ■ d N+1 + • • • + (c a »+«-i), can be written as 
P(di , • • • , d N ) -d N+1 + Q(dt , ••• ,d N ) where P and Q are both poly- 
nomials of the above type in N variables. Then, assuming that the state- 
ment is true for n = N, P + Q ^ and P ■ Q ^ for all positive values of 
the variables d r , • ■ • , d N if and only if at least one of the coefficients 
c , ■ • • , c 2 ,v+,_i is nonzero and all nonzero coefficients have the same sign. 
But, we know that P-d N+ , + Q 7^ for all d N+1 > if and only if P + 
Q^OandPQ^O. □ 

A. I Property (i) is Equivalent to (ii) 

Let D = diag (d, , • • • , d„). By expanding det (AD + B) along the 
first column we have 

det (AD + B) = d x -det P + det Q, 

where the first columns of P and Q satisfy Px = A y , Q^ = Bi , and for 
k = 2, • • ■ , n, P k = Q k = (AD + B) k . Both P and Q are independent 
of (/, . We now expand det P and det Q along their second columns, 
resulting in 

detP = d 2 -det R + det S, 

detQ = d 2 -det U + det V, 

and hence, 
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det (AD + B) = d,d 2 -det R + rf, -det S + d 2 -det U + det 7, 

where 

B L = A, , R 2 = A 2 , 

S 1 = A, , S 2 = B 2 , 

U, = B X , U 2 = A 2 , 

F, = B, , V 2 = B 2 , 

and for k = 3, • ■ • , n, 

R k = S k = U k = V k = (AD + B) k . 

Proceeding in this manner until all columns of (AD + B) have been 
encountered, we obtain an expansion of det (AD + -#) as a polynomial 
in the variables \d t , d 2 , • • • ,d n \ whose coefficients are the determinants 
of the matrices in Q(A, B). By using Lemma 1 it thus follows that (i) 
and (ii) are equivalent. 

A. 2 Property (vi) Follows from (i) and (ii) 

According to (ii) there exists a complementary pair of matrices (M, N) 
taken from Q(A, B) such that det M ^ 0. Let D = diag (d x , • • • , d n ) > 
0, then det (M~ l N + D) ^ if and only if det (MD + N) ^ 0. But, 
using property (i), det (MD + N) = det (AD -f B) ■ det D ^ 0, where 
the matrices I> = diag (d lf ■ • ■ , d„) > and D = diag (d t , ■ ■ ■ , d n ) > 
are defined by d k = d k and 5 fc = 1 if M k = A k , and d k = \/d k , d k = d k 
otherwise (for k = 1, ■ ■ ■ , n). Thus, M~ l N G -Po . 

A .3 Property (i) Follows from (vi) 

Using the notation above, it is clear that for each diagonal matrix 
D > 0, det (AD + B) = det (i¥£ + N) -det 5. Thus, if M"W G Po it 
follows that det (AD + B) ?* 0. 

A .4 Property (vii) is Equivalent to (vi) 

Clearly property (vi) follows from property (vii). Thus, we need 
only prove that (vi) implies (vii). Let (M, N) and (P, Q) both be comple- 
mentary pairs taken from e(.4, B) with M -1 N G -^o and det P ^ 0. 
For any D = diag (d Y , • • • , d n ) > 0, det (P _1 Q + D) ^ if and only if 
det (PD + Q) ^ 0. But det (PZ) + Q) = det (il/Z) + N) ■ det -S ^ , where 
the matrices i) = diag (d x , ■ ■ ■ , d n ) > and D = diag (3i , • • • , d n ) > 
are defined by d k = d k and d lc = 1 if P k = M k , and d k = \/d k , ct k = d k 
otherwise (for k = 1, • • -, n). Thus, P _1 Q G -Po • □ 
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